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ABSTRACT 

The recently discovered system J0651 is the tightest known detached white dwarf (WD) binary. 
Since it has not yet initiated Roche-lobe overflow, it provides a relatively clean environment for testing 
our understanding of tidal interactions. I investigate the tidal heating of each WD, parameterized in 
terms of its tidal Q parameter. Assuming that the heating can be radiated efficiently, the current 
luminosities are consistent with Qi w 7 x 10 10 and Q 2 ~ 2 x 10 , for the He and C/O WDs, respectively. 
Conversely, if the observed luminosities are merely from the cooling of the WDs, these estimated values 
of Q represent upper limits. A large Qi for the He WD means its spin velocity will be slower than 
that expected if it was tidally locked, which, since the binary is eclipsing, may be measurable via the 
Rossiter-McLaughlin effect. After one year, gravitational wave emission shifts the time of eclipses by 
5.5 s, but tidal interactions cause the orbit to shrink more rapidly, changing the time by up to an 
additional 0.3 s after a year. Future eclipse timing measurements may therefore infer the degree of 
tidal locking. 

Subject headings: binaries: close Gravitational waves stars: individual: SDSS 

J065133. 33+284423.3 — white dwarfs 



1. INTRODUCTION 



Compact white dwarf (WD) binaries are important 
as progenitors for a variety of interesting astrophysi- 
cal systems and/or events. Although the merger is 
poorl y understood, i t may resul t in a helium -rich sdB 
star (iSaio fc Jeffervl l200fl 12001 IHeberl l2009h . R CrB 
star (jWebbinkl 11981. survive as mass-tran sferring AM 



CVn binaries (IPostnov fc Y ungclson 2006) which may 
produce a .la supernova (|Bildsten et alj|2007f). or even 
result in a class of Type la supernovae ( Iben fc Tutukovl 
1984). Prior to Roche- lobe overflow, these are among 
the strongest gra vitational wave sources in our Galaxy 
(Nelemans 2009). Tidal interactions are important 
at these short distances, which may alter the inspiral 
and determine the fate upon merger. Unfortunately, 
most tight binar i es are also accret ing (as studied in 
iMarsh et al.ll2004lDelove et alj2005f ). which complicates 
our ability to isolate the effects of tides. 

The detached WD binary SPSS J0 65133.33+284423.3 
(hereafter J0651: [Brown et "all 1201 il l has an orbital pe- 
riod P = 765 sec, the smallest of any such binary yet 
discovered. Furthermore, the system's lightcurve shows 
ellipsoidal variations, Doppler boosting, and primary and 
secondary eclipses, which allow the properties of this sys- 
tem to be tightly constrained. Since there is no accretion, 
J0651 is ideal for studying the role of tides. 

In the following study I explore the effects of tidal in- 
teractions in J0651, focusing on tidal heating and eclipse 
measurements. In Sj2] I summarize the tidal and orbital 
evolution equations and estimate the asynchronicity and 
tidal heating expected. In <J3]I follow the orbital evolu- 
tion with time. I summarize how tides cause deviations in 
the period derivative and eclipse timing from what is ex- 
pected if the evolution was driven merely by gravitational 
waves. This is in stark contrast to the Hulse- Taylor pul- 
sar, which has a period derivative exactly equal to that 
predicted from Einstein's theory ([Hulse fc Taylorlll975l : 



IWeisbereet aDl20loft . In g]I summarize the results and 
discuss future work that can be done on this problem. 

2. GOVERNING EQUATIONS AND ANALYTIC ESTIMATES 

I begin by summarizing the equations used to study the 
time evolution of th e binary. Following the convention in 
iBrown et all ([201 If ), I refer to the less massive He WD 
as the primary, with Mi = 0.25M Q and R\ = 0.035 i? , 
and the more massive C/O WD as the secondary, with 
M 2 = 0.55 M Q and R 2 = 0.013 R Q . 

The tidal force raises a tide on the primary, Hi , induc- 
ing a quadrupole moment 



grav,l 



Hi 
Ri 



AhRj ~ M 2 R\ 



Ri 



(1) 



where a is the orbital separation^. The component of 
the quadru pole that contributes to secular changes in 
the orbit is (|Goldreich fc Soterlll966l ). 
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where t\ ~ l/Qi is the angle of lag due to internal fric- 
tion, ai = 2(n— Qi) is the tidal forcing frequency in the 
rotating frame of the primary (the factor of 2 is due to 
the m = 2 perturbation represented by the tidal defor- 
mation), and n — 2n/P is the orbital frequency. 

The tidal Q that is appropriate for WD binaries is 
not currently known. In other astrophysical systems, Q 
can vary quite dramatically from ~ 10 for rocky bodies 
like the Earth (IGoldreich fc SoterilT96(l to - 10 5 - 10 6 
for gaseous planets such as Jupiter or extrasolar Jupiter s 
(IPeale fc Greenbergll980tlOgilvie fc Lir][200llWull2005[l . 

1 For these derivations I focus on the primary, but these equa- 
tions can just as well apply to the secondary by switching the 
subscripts 1 and 2. 
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For WD binaries, [Campbell (|1984l ) estimated synchro- 
nizati on times of ~ 10 8 yr , which correspond to Q ~ 
10 12 . IWillems et all ([27)10 ) found Q ~ 10 15 , but their 
analysis did not include dynamical tide effec ts, so this 
is merely an upper limit. iFuller fc Lail (120111) looke d at 
dynamical tides and then, appealing to iWul (|1998[ ) for 
damping times, estimate Q ~ 10 4 — 10 11 (depending the 
radial waveleng th of the dom i nant excited mode). Other 
studies, such as lRacine et al.l (|2007l) . are more focused on 
when accretion is occurring. Due to these uncertainties, 
I take Q to be a free parameter for the present study. 
This secular quadrupole moment implies a torque 
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Assuming solid-body rotation, the spin evolves as 



dt 



(iifii) ~ Ax 
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(4) 



where h = kM\R\ and I set k = 0.2 (jMarsh et al.ll2004 
but this may be modified for the hot He WD). 
The angular momentum of the binary evolves as 



J, 
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(5) 



w here the angular momentu m loss to gravitational waves 
is ([Landau fc Lifshitdl 19751) 
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where J orb = (Ga/M) 1 ^ 2 MiM 2 and M = M x ^ 
The evolution depends on two key timescales. 
is the gravitational wave timescale, 
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dP 



dt 

= 5.8 x 
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where M . 25 = Mi/0.25M Q , M . 55 = M 2 /O.55M , 
Mo.s = M/0.8M Q and P 3 = P/10 3 s. The other is the 
torquing timescale found from equation (U), 



n = ai 



= 1.6 x lO^o^gioMo^Mo.ssMo.g^^P^ yr, (8) 

where P9.3 = R\/2 x 10 9 cm, fc . 2 = fc/0.2 , and Qiq = 
Qi/10 10 . In previous calculations, such as IMarsh et al.l 
(2004) , the strength of the tidal torquing is discussed in 
terms of a tidal synchronization timescale, r s . This is 
merely a factor of 2 smaller than the T\ defined here. 
This demonstrates that if Qi is assumed constant, the 
synchronization timescale decreases rapidly as the binary 
shrinks. 

Given these timescales, I estimate the asynchronicity 
of the primary's spin. The time derivative of the tidal 
forcing frequency is 



da 1 



dt \ dt 



~~dT 



= 2 



Tl 



(9) 



The competing effects of the gravitational wave emission 
promoting asynchronicity and the torques tidally locking 



the stars drive the system toward an equilibrium where 
da±/dt w 0. The tidal forcing frequency is then 
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The general trend is that as the He WD gets closer to 
the C/O WD, the tidal torques become stronger (since 
Ti oc P 3 and r gw oc P 8 / 3 ), which in turn makes ai/n 
smaller. On the other hand, for a larger Qi there is 
less of a lever arm for torquing the He WD, and its spin 
becomes more asynchronous. 

A tidal torque acting on an asynchronously spinning 
WD implies that work is being done. The rate of energy 
input is 
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Qi 



= 2.3 x 10 31 fc 2 2 ( 
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(11) 



Whether or not this E\ is observed depends on how 
shallow it is deposited. One expects that the luminos- 
ity of the He WD is roughly given by E\ when the 
thermal timescale down to where the tide is damped 
is less than r gw . Ilben et aT] (|1998[ ) consider this issue 
in more detail by distributing the tidal heating evenly 
thro ughout the WD, an d then following the WD's cool- 
ing. IBrown et al.l (|2011[ ) estimate the effective tempera- 
tures for the primary and secondary WDs to be 16, 400 K 
and 9, 000 K, respectively, providing luminosities of L\ — 
3.1 x 10 32 ergs s" 1 and L 2 = 3.8 x 10 30 ergs s^ 1 . Using 
equation (fTTj) . Q\ w 5 x 10 10 is needed for E\ w L\. 
Equation (jlip appears to imply that the heating in- 
creases indefinitely with larger Qii but the equilibrium 
condition only applies for t\ < r gw . Using a maximum 
asynchronicity frequency of o\ = 2n, I find 

Pl.max » 5.9 X 10 32 Qr 1 ^0.25M "8 1 ^9 3 P 3" 5 er S S 

(12) 

for the maximum heating rate. 

The above estimates focus on the asynchronicity of the 
He WD, but the same arguments can be extended to the 
C/O WD. Using equation (Tig)) . 
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where the term on the far right assumes the masses and 
radii appropriate for J0651. This result reflects the fact 
that the tidal torque is smaller on the C/O WD, so the 
asynchronicity is larger for the same Q. The ratio of the 
heating rates is 



P2 



Q2 (M2 
Qi \ Mi 



Pi 
P 2 



29-^. 



(14) 



If the observed luminosity of the C/O WD is also due 
to tidal heating, then Q 2 w (Q l /29)(L 2 /L 1 ) « 2 x 10 7 . 
The structure of the He and C/O WDs are quite differ- 
ent, with the former having a shallower convective region. 
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In future work it would be interesting to theoretically ex- 
plore whether such differences lead to Q2 "C Qi- 

Numerical simulations of merging WDs often assume 
that t he WDs are not tida l ly locked at merger (for ex- 
ample [SegretainetaDIillS [GuerreroeraD[200l) . From 
our derivation we can see that this is not always the case, 
and in general it depends on the value of Q. For example, 
^l/n ^ 0.1 merely require Q x < 5 x 10 9 . If Q 2 = 2 x 10 7 
right up until merger, then 02/71 « 0.02 and the C/O 
WD must be tidally locked unless Q2 increases by a cou- 
ple of orders of magnitude as the binary nears merger. 

In exploring th e survival of binary WDs in forming 
AM CVn systems, iMarsh et all (|20o4 conclude that the 
synchronization time for the secondary must be less than 
1, 000 yr. For the model presented here, that would imply 
Qi ^ 10 s when accretion first begins. The current L2 
appears to place a more stringent limit on Q2 than this 
already, although this conclusion again depends on how 
Q2 changes as the orbit shrinks. 

3. NUMERICAL INTEGRATIONS 

In Figure Q] I show numerical integrations of equa- 
tions (jlj and (0, using the mass and radii appropriate 
for J0651. The tidal Q parameters are set to be constant 
at Qi = 7 x 10 10 and Q2 = 2 x 10 7 , so as to give heat- 
ing rates at P = 765 s that are the same as the present 
luminosities of each star. The WDs are assumed to be 
non-spinning initially at a large orbital period, but are 
found to be quickly spun up by tides until da/dt « is 
reached, consistent with the assumptions for my analytic 
estimates. The integration ends when the Roche-lobe 
around the He WD becomes equal to its radius, which 
occurs at P » 420 s (ignoring potential changes to i?i 
due to tidal heating). The He WD is spinning signifi- 
cantly more slowly than the orbital period because of its 
large Q\. The vertical dotted line denotes the current 
location of J0651 at 800, 000 yr before merger. If the 
values of Q remain constant, the luminosity of the He 
WD increases by a factor of ~ 15 before tidal disruption. 

In the top panel of Figure [2 I calculate the rotational 
velocity of the primary V% = as a function of or- 

bital period. When the values of Q are chosen to match 
the current luminosities, V\ ~ 120 km s _1 . Another case 
where Qi = Q2 = 10 7 is also plotted, which represents 
what happens when the WDs are nearly tidally locked, 
resulting in V\ ~ 200 km s . Since the binary is eclips- 
ing, the difference between these cases may be meas ur- 
able via the Rossiter-McLaughlin effect (|GrootJl201rr) . 

An important consequence of the tidal interactions is 
that the orbital period derivative deviates from what is 
expected if the system is purely driven by gravitational 
wave losses. Using equation (JT]), gravitational waves 
alone give 



ip , 



1.7 x 10- 4 A/ .5 5 M . 25 M j /3 P, 5/3 



s yr 



(15) 

so for J0651 one would expect a period decrease of 
|P gw | = 2.7 x 10~ 4 s y^ 1 . When tidal effects are in- 
cluded, the period decreases faster, because angular mo- 
mentum is loss by the orbit to spin up the WDs. I assess 
this effect in the numerical models using the relation 
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Fig. 1. — The binary evolution as a function of time, using Qi = 
7 X 10 10 and Q2 = 2 X 10 7 . Masses and radii are chosen to match 
J0651. The top panel shows the spin period of each star. The 
orbital period of the binary is nearly equal to the spin period of the 
C/O WD (dashed line) and thus is not plotted. The middle panel 
plots the tidal heating rate for each star, and the bottom panel 
shows the surface effective temperatures. The vertical dotted line 
shows the current location of J0651. 




» 6 



0.1 



■a±. 0.05 













1 


II 

Q 1 = 7xl0 


,0 , Q 2 = 2xl0 7 



Only GWs 




900 800 700 600 500 400 

P (s) 

Fig. 2. — The top panel plots the velocity of the primary 
Vi = Q1R1 as a function of orbital period. I consider a case where 
the values of Q are chosen to fit the observed luminosities (solid 
line) , and a low Q example where the WDs are nearly tidally locked 
(dashed line). The middle panel plots the period derivative, com- 
paring merely gravitational wave losses (dot- dashed line) with the 
other two tidal cases. The bottom panel shows the deviation of the 
two tidal cases from pure gravitational wave losses. The vertical 
dotted line shows the current period of J0651. 
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(16) where J orD is found from equation (J5j). 
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In the middle panel of Figure [H I compare the pe- 
riod derivative for purely gravitational wave losses (dot- 
dashed line) with two different degrees of tidally lock- 
ing (solid line and dashed line). This confirms that the 
period decreases more rapidly when tidal effects are in- 
cluded, and this deviation increases as the binary inspi- 
rals to smaller orbital periods. 

In the bottom panel of Figure [2 I plot the fractional 
change in the period derivative, comparing the tidal cases 
with purely gravitational wave losses. For the values of Q 
that best fit the current luminosities, the period deriva- 
tive is » 3.0% larger, whereas if the WDs are nearly 
tidally locked, the period derivative changes by « 5.6%. 
Looking for these deviations in future observations would 
provide an important measurement of the presence and 
level of tidal locking. An observed \P\ > \P gw \ implies 
strong tidal locking, but if \P\ sa |-Fgw|, then the asyn- 
chronicity and tidal heating must be large. 

A non-zero P can also be measured by its affect on the 
eclipse timing. After a time t, this changes by 

t 2 P 

Change in eclipse time = (17) 

For purely gravitational wave emission, the change in the 
eclipse time is 5.5 s after one year. When tidal interac- 
tions are included, this change is 3—6% larger, which cor- 
responds to eclipses 5.6 — 5.8 s sooner. Associated phase 
shifts would be important in studies of this binary by 
the proposed Laser Interferometer Space Antenn a mis- 
sion (|Webbink fc HanlH99l iNelemans et al.ll2004[ ). 

4. DISCUSSION AND CONCLUSION 

I considered the effect of tides in the detached binary 
J065I. Assuming that the current luminosity of each 
WD reflects its tidal heating, I found Qi = 7 x I0 10 
and Q 2 = 2 x I0 7 for the He and C/O WDs, respec- 
tively. The values of Q cannot be greater than this, oth- 
erwise the tidal heating would be inconsistent with the 
current luminosities. The degree of tidal locking can be 
measured in future observations, either via the Rossitcr- 
McLaughlin effect or eclipse timing. If only gravitational 
waves are acting, then after one year the eclipses occur 
5.5 s earlier. Tidal locking causes this difference to in- 
stead be as much as 5.8 s. Since the time of eclipses 
scales oc t 2 (see eq. [H]), longer baseline observations 
will be important for making an accurate measurement. 

Before deviations of P from merely gravitational wave 
losses can be inferred, the masses of the WDs must 
be measured sufficiently accurately. Since P gw oc 
M1JW2/M 1 / 3 (eq. [15]). an uncharacteristically large \P\ 
may instead imply larger masses. Mass measurements 
within a couple of percent would be ideal, which in turn 



means that the inclination must be highly constrained. 
At the same time, the P cannot be strictly due to gravi- 
tational waves, otherwise the tidal heating would be too 
large in comparison to the observed luminosities. Once 
P is measured in future observations, simultaneous mod- 
eling of the masses and luminosities should be done to 
constrain what values of Q, and thus degree of tidal lock- 
ing, are allowed. 

A critical question about J0651 is its age and how it 
relates to the larger population of WD binaries. The 
current merger time for J065I is 800, 000 yr, while the 
age of the He WD inferred from its current effective 
temperature pres uming no tidal heating, is ~ 10 8 yr 
(|Panei et al.l 120071 ). On the face of it, this would im- 
ply that there should be « 100 binaries with periods 
of P(r gw = I0 8 yr) w 50 min for every system like 
J065I . Although we are limited by small number statis- 
tics, such a large number of binaries has not been found 
(jKilic et al.H20lTh . This may point to an incompleteness 
of the current surveys. On the other hand, the progeni- 
tors to J065I may be wider binaries that are older, cooler, 
and harder to detect. As the binary contracts, the tidal 
heating then makes the He WD bright and easier to ob- 
serve. A more self-consistent model of the tidal heating 
plus stellar cooling, including changes to the He WD ra- 
dius, is needed to assess what luminosity is expected as 
a function of time. Also, depending on when the heating 
takes place it can either be radiated readily or have time 
to heat the core and alter the stellar structure. I plan to 
explore these details better in a subsequent study. 

Although there has been some work on estimating the 
Q of WD binaries (see iFuller fc Laill201ll and references 
therein), more can be done, especially in the context 
of He WDs. The presence of a thin layer of hydro- 
gen (~ I0 -3 Mq) on the He WD will alter the eigen- 
functions of the excited modes, determining which spe- 
cific modes are driven as well as affecting their thermal 
damping timescale, which in turn determines Q. In a ddi- 
tion, the lightcurves presented bv lBrown et al.l (|2011l) for 
J065I show considerable dispersion (see the upper panel 
of their Figure 4) that could potentially point to stellar 
oscillations. The hot temperature of the He WD makes it 
well outside of the traditional instability stri p associated 
with WDs, which requires T eff < 11,000 K (|Arras et al.l 
120061 ; fStcinfadt et al. 2010). But strong tidal interaction 
may be another method of exciting observable oscilla- 
tions, which should also be investigated. 
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grant NNX07AH06G, NSF grant AST-0855535, and by 
the Sherman Fairchild Foundation. 
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